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Abstract 

We have conducted multi-epoch synthesis imaging of A =2 and 3 millimeter (mm) continuum emission 
and near infrared K band (2.2 /xm) imaging of a flare event in January 2003 that occurred on the young 
stellar object GMR-A which is suggested to be a weak-line T Tauri star in the Orion cluster. Our mm data 
showed that the flare activity lasted at least over 13 days, whereas the K-band magnitude did not change 
during this event. In addition, we have succeeded in detecting short time variations of flux on the time 
scales of 15 minutes. The total energy of the flare is estimated to be ~ 10 35 ~ 36 erg, which makes it one 
of the most energetic flares reported to date. Comparing the mm continuum luminosities with reported 
X-ray luminosities, we conclude that the mm flare was similar in nature to solar and other stellar flares. 
Our results will be a crucial step toward understanding magnetically induced stellar surface activities in T 
Tauri stars. 
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1. Introduction 

Recent radio and X-ray observations have revealed that 
pre-main-sequence stars show flare activity as seen in the 
Sun and late-type stars. The radio emission has been 
attributed to gyrosynchrotron emission from mildly rela- 
tivistic electrons and is preferentially seen in X-ray bright 
stars. Luminosities of these flares at centimeter (cm) 
and those in soft X-rays are known to correlate well with 
each other (Giidel 2002; references therein). Moreover, 
Feldman, Laming & Doschek (1995) discovered a correla- 
tion between the peak temperature of a flare (T) and its 
volume emission measure [EM) derived from X-ray spec- 
tra of solar and stellar flares. Shibata & Yokoyama (1999, 

2002 hereafter SY02) attempted to explain the EM - T 
relation by magnetohydrodynamic models. 

Stellar X-ray emission is highly time variable, but vari- 
ability at millimeter (mm) wavelengths is poorly un- 
derstood. This Letter 1 , presents in detail our mm in- 
terferomctric observations of a flare that occurred in 

2003 January in GMR-A. Initial results were reported 
in Nakanishi ct al. (2003). The flare was originally 
discovered by Bower, Plambeck & Bolatto (2003a) us- 
ing the BIMA array on January 20 UT at 86 GHz and 
subsequently reported by Bower et al. (2003b, hereafter 
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BPB03). GMR-A is identified as a member of the Orion 
cluster (M42; distance of 450 pc) from early VLA ob- 
servations (Garay, Moran & Reid 1987). BPB03 identi- 
fied that the object is a K5 type weak-line T Tauri star 
(WTTS; bolometric luminosity ~ 6L©) by near-infrared 
(NIR) spectroscopy. During the flux monitoring campaign 
at cm wavelengths over 7 months in 1990, GMR-A showed 
one radio flare (Felli et al. 1993). 

2. Observations 

Synthesis imaging of the A —2 and 3 mm continuum 
emission was carried out using the Nobeyama Millimeter 
Array (NMA) with the most expanded configuration 
(Tabic 1). Note that the two frequency bands were not 
observed simultaneously. The phase tracking center was 
set at R.A.=5 h 35 m ll?8, Decl.=-5°21'49"2 (J2000) and 
the field of view (FOV) was 70" at 3 mm and 50" at 2 
mm. For the backend, the 1 GHz bandwidth correlator 
was employed. We used the continuum source 0420—014 
as a bandpass calibrator and 0528+134 as a phase and 
gain calibrator. From observations of Uranus, we esti- 
mate the fluxes of 0528+134 to be 2.56+0.5 Jy at 98 GHz 
and 2.11+0.4 Jy at 147 GHz. 

The NIR K-band ( A =2.2 /jm) data were taken with the 
Okayama Astrophysical System for Infrared imaging and 
Spectroscopy (OASIS) mounted on the 1.88 m telescope 
of the Okayama Astrophysical Observatory OASIS has a 
FOV of 4' with an image pixel size of 0"97. The obser- 
vations were done with seeing of J$l."7 on 2003 January 
18.451 and 21.459 UT and its calibration was performed 
using 23 stars reported in Muench et al. (2002). 



2 R. S. Furuya et al. [Vol. . 



Table 1. Summary of the NMA Observations toward the GMR-A Radio Flare in 2003 January 



UT 2003 


Time Ranj 


3 c of HA* 


Timc^ 


Frequency 


Sensitivity-! 1 


n v a § 
^maj * "min 




(hh:mm:ss.s) 


(hrs) 


(GHz) 


(mjy beam -1 ) 


(arcsec x arcsec) 


Jan.24.56 


02:02:29.9 


03:15:37.9 


1.00 


97.782 


4.4 


2.51x1.18 


Jan. 25. 55 


00:31:08.8 


01:45:22.1 


1.00 


146.969 


5.1 


1.32x0.91 


Jan.25.61 


02:30:29.5 


03:18:37.5 


0.750 


97.782 


5.4 


2.28x1.21 


Jan.28.44 


-02:46:34.1 


-00:46:13.7 


1.90 


98.243 


4.1 


2.25x1.35 


Jan.28.58 


00:42:59.0 


03:02:22.6 


1.75 


141.359 


2.4 


1.65x0.92 


Feb.02.57 


00:32:38.8 


03:22:05.5 


2.25 


99.023 


4.1 


2.89x2.50 


Feb.07.57 


01:08:22.9 


03:41:48.3 


2.00 


99.023 


2.8 


2.58x1.35 


Feb. 15.37 


-04:02:00.0 


00:07:37.6 


3.25 


99.023 


3.1 


2.85x1.56 



* Time ranges of the observations given with starting (left) and ending (right) hour angles (HAs), ' Effective integration 
times, * RMS image noise levels, § Synthesized beam sizes along major (0 ma j) and minor (# m i n ) axes. 
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Fig 

map (contour) taken on 2003 January 25 and the OASIS 2.2 
fim image on January 18. Solid and dashed contours indicate, 
respectively, positive and negative ones starting at 2tr level 
with ±ltr step of the noise level (Table 1). The cross indicates 
the peak position of the 86 GHz continuum emission (BPB03). 
The synthesized beam size of the NMA observations is shown 
in the bottom right corner. 

3. Results and Discussion 

3.1. Overall Properties: The Maps and Time Variations 

Fig. 1 is a map of the continuum emission at 2 mm 
superposed on a 2.2 /urn image. The peak positions 
of the 2 mm and 2.2 /jm emission coincide with each 
other within errors (< 1"), and their positions agree with 
the previous measurements with BIMA and VLA (e.g., 
Felli et al 1993). The K band flux was measured to 
be 9.625±0.032 mag on January 18.451 and 9.586±0.047 
mag on January 21.46. There was no significant flux 
change at K band during the most prominent mm flare on 
January 20.188 (BPB03). Furthermore, the above mea- 
surements arc consistent with those in 1997 (Muench et 
al. 2002) and in 1999 (Hillenbrand & Carpenter 2000) 
within 0.05 mag, suggesting that the K band flux is sta- 
ble. Comparing our K band data with those of H band 
(Hillenbrand & Carpenter 2000), GMR-A is embedded in 
dust {{H - K) = 2.311±0.047 mag). Such strong infrared 
excess is more consistent with a classical T Tauri star 
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Fig. 2. Time variations of (a) the mm continuum emission 
fluxes and (b) its spectral indices between 2 and 3 mm. In the 
panel (a), open and filled circles indicate 2 and 3 mm fluxes, 
respectively. The horizontal dashed lines in (b) indicate the 
possible range of free- free emission (Rodriguez et al. 1993). 

rather than a WTTS, as BPB03 argued as well. 

Table 2 summarizes flux densities measured at 2 and 
3 mm integrated over one day. We present radio "light 
curves" in Fig. 2a where a clear pattern of a flare can be 
seen at 3 mm: the flux density at 3 mm increased from 
January 25 to 28, then it gradually decayed and was not 
detected on February 15. Although BPB03 reported that 
the onset of the burst was January 20, we did not de- 
tect any emission at 98 GHz on January 24. BPB03 also 
reported non-detections in their VLBA observations on 
January 24 at 15 and 22 GHz. We suggest that these non- 
detections could be due to intrinsic short time variations 
of the source (§3.2) or a possibility that the January 25 
detection is attributed to a different (possibly the second 
giant) burst event as BPB03 suggested. It is also possible 
that January 28, February 2, and 7 measurements refer 
to independent flares that may occur frequently on the 
star. The behavior at 2 mm, on the other hand, was sig- 
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Table 2. Summary of the Millimeter Continuum Emission Measurements 





3 mm 


2 mm 


Spectral 


UT 2003 


Frequency 


Flux 


Frequency 


Flux 


Index 




(GHz) 


(mJy) 


(GHz) 


(mJy) 




Jan.24 


97.782 


< 11.3* 








Jan. 25 


97.782 


22.3±5.4 


146.969 


39.1±5.1 


+ 1.39 ±0.78 


Jan.28 


98.243 


36.2±4.1 


141.135 


< 6.3* 


< -2.4t 


Feb.02 


99.023 


23.3±3.5 








Feb.07 


99.023 


18.2±2.8 








Feb. 15 


99.023 


< 8.0* 









* 3cr upper limit, t 5<r upper limit. 
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Fig. 3. A series of 2 mm continuum emission maps taken in 
2003 January 25 with a 15 min bin. The UT shown in each 
panel is the center time of each integration. The ellipses in the 
bottom right corners indicate the NMA beam sizes. Contour 
intervals are the same as in Fig.l and the typical noise level 
is 12 mJy beam -1 . 

nificantly different: though it showed bright emission on 
January 25.55, it was not detected at January 28.58. 

The above time variation leads to a drastic change in the 
apparent spectral indices. We derived indices between 2 
and 3 mm assuming a power-law type spectrum (S v oc v a ): 
the index, a, observed on January 25 was positive with 
+ 1.39 ±0.78, but the January 28 index was negative with 
a 5cr upper limit of —2.4. The former could be explained 
with a combination of optically thin free-free emission and 
thermal emission emanating from dust, but the latter can- 
not be easily explained. Rodriguez et al. (1993) reported 
that a negative spectral index less than —0.1 cannot oc- 
cur for thermal processes involving free-free emission and 
absorption (Fig. 2b), regardless of the electron density and 
temperature distribution. In fact, we found that the nei- 
ther of the indices can be reproduced with the model of an 
H II region embedded in a dust cloud core by Furuya et al. 
(2002). Cool dust radiating at those wavelengths should 
not exhibit such rapid variations of a, as demonstrated by 
our K band observations. 
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Fig. 4. Plots of mm continuum fluxes vs. PAs (see §3.2) 
at (a) 2 mm and (b) 3 mm. Each point is obtained with 
15 min integration with an interval of 20 min. (a) The 2 
mm data were taken on January 25: small labels show the 
corresponding panel numbers in Fig. 3. (b) The 3 mm data 
shown with open circles and triangles are taken on January 
28 and February 2, respectively. The dotted horizontal lines 
show typical image noise levels. 



3.2. Short Time (< A Day) Variations 

Solar and stellar flares are commonly accompanied by 
short time scale (^ 1 day) variations at cm radio and X- 
ray wavelengths (Dulk 1985). In fact, BPB03 detected 
short time variations at 86 GHz within a day of January 
20. To investigate such rapid variation, we selected data 
sets containing 15-minute (min) integrations in which the 
emission was detected with signal to noise ratios higher 
than 3: the 2 mm data taken on January 25, and the 
3 mm data on January 28 and February 2. Fig. 3 is a 
series of 2 mm continuum maps taken every 20 min on 
January 25. The highly time variable nature of the source 
can clearly be seen. At the beginning, its flux was 46±11 
mJy (Fig. 3a), it flared up to 68 ± 10 mJy after 20 min 
(Fig. 3b), then it was decaying (40±9 mJy; Fig. 3c) and di- 
minished (3cr upper limit of 33 mJy: Fig. 3d) 60 min from 
the beginning. During the observations, the calibrator 
was observed every 20 min with 4 min integrations, and 
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the amplitude of its raw visibilities taken by each baseline 
during each 4 min showed no systematic variation with 
fluctuation less than 12% of the mean value. Therefore 
the observed flux variations are not an artifact of calibra- 
tion. 

It is important to discuss how linearly polarized emis- 
sion affects the flux measurements because radio emission 
showing short time variations might be intrinsically lin- 
early polarized. Receiver systems detect one direction of 
incident radiation. The NMA receivers on each antenna 
are placed at the Nasmyth foci for each band and are 
aligned so that detected directions with individual anten- 
nas are parallel when they are projected to the plane of 
sky. Since the direction rotates while tracking the ob- 
ject, observing with such systems could cause apparent 
flux variation. Assuming that the radiation has a lin- 
early polarized component and that its polarization angle 
is constant, the phase of the flux variation should be syn- 
chronized with the position angle (PA) of the receiver with 
a 180° period. To verify lack or existence of such a peri- 
odic pattern, we present Fig. 4 where the flux variation is 
transformed function of PA. 

The flux variations of the three data sets in Fig. 4 do 
not show any periodic pattern of phases. The 2 mm fluxes 
(Fig. 4a) show variation on a time scale much shorter than 
a 180° period. This strongly suggests the 2 mm flux vari- 
ation is real. The 3 mm data taken on January 28 seem to 
show a considerable flux change peaked around PA=— 30° 
but there is no signature which suggests an existence of a 
180° period (Fig. 4b). The other 3 mm data contain too 
many non-detection points to judge whether there is a si- 
nusoidal variation with a 180° period. Moreover, we did 
not find a unique solution that all the 3 mm data points 
satisfy with a 180° period, even if we consider the overall 
decay of the flux (Fig. 2a). We thus conclude that there is 
no significant contribution of linear polarization and that 
the observed short time variation with time scale less than 
an hour is due to intrinsic variations of the source. These 
results are consistent with those from the BIMA measure- 
ments. 

3.3. The Nature of the GMR-A Flare in 2003 January 

Fig. 5 compares radio and soft X-ray luminosities of 
GMR-A in quiescent and flaring phases with those of 
magnetically active stars in quiescent phase (Giidel 2002). 
Solar flares which are less energetic than the stellar flares 
presented here are known to satisfy the correlation as well. 
Considering that GMR-A is most likely to be a single 
WTTS (BPB03), we conclude that the flux enhancements 
in 2003 January were due to stellar flares from a magnet- 
ically active WTTS, such as V773 Tau (e.g., Feigelson & 
Montmerle 1985). This is because the locus of GMR-A 
in quiescent phase is basically consistent with the correla- 
tion. 

Physical parameters for a flare can be estimated from 
the data of January 20 on which the radio and X-ray ob- 
servations were carried out simultaneously. Using the soft 
X-ray luminosities in the flaring phase and the time scale 
of the mm flare (BPB03; also our results from §3.2) to- 
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Fig. 5. Comparison of Lx and Lr of GMR-A in its quiescent 
and flaring phases with those for magnetically active stars 
in quiescent phase (Giidel 2002). The luminosity ranges for 
GMR-A were calculated using the quiescent phase cm radio 
fluxes reported in Felli et al. (1993), quiescent and flaring 
phase soft X-ray fluxes in Getman et al. (2003) and the flaring 
phase radio fluxes shown in Fig. 4. 

gether with the reported time scale of the X-ray flare (a 
few 0.1 day; Getman et al. 2003), we estimate total en- 
ergy (E) of the flare of ~ 1 35 - 36 erg, taking L x of 10 32 
erg s" 1 (Getman et al. 2003). Once we know E, we can 
pinpoint the locus of the flare in the EM — T diagram (see 
Fig. 8 in [SY02]): we estimate T of ~ 5 x 10 7 K and EM 
of ~ 10 53 cm~ 3 . Subsequently, using Fig. 2 in SY02, we es- 
timate a magnetic field strength (B) of 50 — 150 G, which 
indicates that the cyclotron frequency should appear in 
the GHz range. This estimate is consistent with measure- 
ments of the stellar photosphcric magnetic field strength 
(a few kG; BPB03) which is expected to be an order of 
magnitude stronger than those of coronal magnetic field 
strengths. 

We believe that the T and B estimated above are appli- 
cable to the flares detected with the NMA. This is because 
the observed flux density (Table 2) and the flare's size 
scale that can be estimated from plausible Alfven velocity 
(0.05c assuming B =100 G and typical solar flare den- 
sity of 10 9 cm~ 3 ) and the flaring time scale give us the 
brightness temperature of ~ 10 8 K, which is consistent 
with the T estimated above. In conclusion, our results 
are an important step toward understanding magnetically 
induced stellar surface activities in T Tauri stars and will 
add to the growing evidence that magnetic rcconncction 
processes could explain solar and stellar flares over various 
energy scales. 
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